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We have veri fied that the variati onal mean field theory approach suggested by Narimanov and 



Varma (preprint cond-mat/0002191) being applied to the realistic two-band model provides a good 
agreement with experimental data for the Curie temperature in doped manganites Ai_i,Ba:Mn03 
{x ~ 0.3). We have also considered the problem of an interplay between the ferromagnetic and 
antiferromagnetic interactions by using the same approach. 
The problem of finding the value of the critical temperature for the paramagnetic-to-ferromagnetic transition in 
doped manganites has presented a lot of interest recently. An agreement of theoretical predictions with experimen- 
tal results for the Curie temperature may shed some light on physical mechanisms driving the transition near the 
concentration region, where the colossal magneto-resistance (CMR) effect takes place. 

An interesting approach for obtaining the value for the transition temperature has been recently proposed by 
Narimanov and Varmao (sea-|alsd3) and is based on the variational mean field theory within the formalia™ of the 
double exchange (DE) modefro. The final result for the critical temperature dependence on concentrationoH can be 
approximately written in the following form: 

T, ^x{l-x) W, (1) 

where W is a bandwidth. For W = 1.8 eV and 0.15 < x < 0.3 Eq. (|l]) recovers the experimentally observed values for 
the Curie temperature. This resuit-,can be generalized in such a way, that the DE model works well in order tCj-describe 
the physics near the CMR regioEO'El and in principle one does not need to take into account other mechanismsB'El which 
might effect the transition. l— . 

In what follows we adopt the variational mean fieid-theory scheme for the DE model via degenerate orbitalsEEl in 
order to verify the validity of the two-band approactO'EJ near the CMR region. In the DE model within the two-band 
scheme, there is no electron-hole symmetry, which means that the Curie temperature is non-symmet |ri j | w ith respect 
to substitution a; — > 1 — x. Since for the small concentrations percolation effects play an essential roleMlij, we will be 
primarily concerned with obtaining the value for the Curie temperature near the CMR region [x ^ 0.3). As we will 
see, the approach we are using works well and recovers the right numerical value for the Curie temperature in units of 
hopping amplitude the value for which in the two-band model for manganites has been estimated inll3. The relative 
simplicity of present approach-alaa allows us to consider a problem of suppression of the ferromagnetic transition by 
the superexchange interactionliJilj. 

We consider three t2g electrons on a Mn ion as the localized classical spins. The Hamiltonian for the DE model on 
the degenerate orbitals can be defined as: 

HdE = - ^tl^-cl^a^jba - JH^S,-c\^^aa(3C,a!3- (2) 
(ij> i 

The matrix elements describe the electrons, hopping from one site to another, on two-fold degenerate Cg orbitals. 
As it was shown by Anderson and Hasegawa @, in the limit of strong Hund's coupling, the hopping amplitude of the 
electrons between the sites i and j acquires an angular dependence ty = t{9ij), where 9ij is the angle between the core 

spins Si and Sj. In the DE model the hopping amplitude is the largest, when the spins are aligned, which means that 
a ferromagnetic state minimizes the kinetic energy. Thus, the effective Hamiltonian for our problem can be written 
as: 

H = Hde + (3) 

We employ the variational mean field theory by introducing the spin orientation distribution function P{^) = 
^P('(?). In the mean field approximation, the distribution function depends only on single spin orientation, which is 
described by an angle The idea is to represent the free energy of our system as a functional of P('i?) and, following 
the standard mean field treatment, to find the self-consistent equation for the distribution function. 

The free energy is defined as: 

-P ^el ^" -^spin P ' Sspin : (4) 
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where Sspm is the entropy of our system, E^i is electron energy, E'spm is the spin energy. 

As mentioned above, the distribution function for the system of spins is just the product of the distribution functions 
of the individual spins. Thus, in the semiclassical limit, the entropy can be written as: 

Ssp^n = - J d^sm^P{^)log[P{§)]. (5) 



In the absence of a magnetic field, the spin energy contribution, Espm, in the classical limit, is just the average of 
the superexchange Hamiltonian: 



Z7T ZTT TT TT 



(6) 



where cos(6') — cosi?i cosi92 + sin?9i sini?2 cos((/)i — (j)2)- 

The electron energy is calculated as follows. First, we calculate the energy for a given spin distribution, E[t] and 
then we just take the average of E[t] over all possible spin configurations: 

E,i = j dtP[t)E[t\. (7) 

To calculate E[t\ we adopt the tight binding approxiKiation. For the two-dimensional representation eg we choose 
the normalized complex functions ipi and given bjiiil: 

ipi<x-? + ex^ + e^y"^, ^2 = V'l*, (8) 

where e = exp(27ri/3). 

The cubic spectrum consists of two branches: 

£i,2(p) = -t{0) • |(c^ + Cy + c^) ± \Jcl + cl + cl- c^Cy - CyCz - c^Cx I (9) 

(we introduced the notations Ci = cos(/eia), i = x,y, z). Since there is an interaction between the Cg electrons and the 
magnetic background, following the discussion by Anderson and Hasegawa, we have introduced the effective hopping 
amplitude t{9) —\ A \ cos{9/2) [9 is the angle between the two neighboring spins and | A | is a hopping amplitude, 
whose numerical value will be defined later). 

Thus, the electron energy for a given spin configuration is: 
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- E / 7^=^(k' 0% (10) 

i—l ^ ' 

where np is Fermi distribution function. 

Using the expressions (|,|,|l^) and the La grange variational method, the free energy functional acquires the following 
form: 

2tt 2-k TT 7r 

i^[i^>5;M,e] ^ J^J^J d^isin^iPi^i) J d^2sin^2P{^2) x 
000 

J^i£^{K0)-X)xnF[e^ik,9)] + ^cos{9)^+T j ddsm^x 

TT 27T TT 

P(t?) log[F(z9)] + C J di^sin^Pi^) ~ fisSB J ^ J ^ 

00 
sini9iCosi?iF(i9i) +-F[A,C] (11) 

(A is a Lagrange multiplier, being a constraint to have a constant number of conduction electrons in our system, ^ 
is a Lagrange multiplier, which provides the normalization for P{'&) and _F[A,^] is a independent part of the free 
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energy, i? is a external magnetic field) . Since the number of conduction electrons is conserved in our system, the free 
energy functional acquires its minimal value for A being equal to the chemical potential /i, which is defined, in our 
case, by the following expression: 



i=l 



(2^ 



nF[{e,ik,e = 0)~^,)/T], 



where x is the hole concentration. 

Now, taking the functional derivative of (O) with respect to f we obtain the following integral equation: 



P(t?) = exp 



E 



(27r)3 



T 



nF[Si{k, 



2T 



cos{0) 



e/T+^cos^' 



(12) 



The integral equation ( |l^ ) is nothing but variational mean- field theory equation for the distribution function P{'d), 
where Pji!}) plays the role of the order parameter. If the system is in the paramagnetic phase, PI'S) is uniform. 
Equation ( |l2|) allows a numerical solution. Its results are shown on Fig. 1. 




FIG. 1. The single spin orientation distribution function P(i9) is shown as a function of i9 for different values of the antifer- 
romagnetic exchange constant J. The results are obtained for hole concentration x ~ 0.3. 

As soon as the temperature decreases (in the presence of a small magnetic field) and approaches the Curie temper- 
ature Tc, small deviations from the miiform distribution will appear in P(i?): 



miT-.T. = o + A(^?), A(^)«l. (13) 



Since A(??) is small, we can expand both sides of Eq. ( [l^ ) in powers of A(i?), keeping only the linear terms and 
solving the final expression for A(i?). The results can be written in the following form: 



. , „, B cos t9 
A(^) oc 

T, = I dd sin I? cost^l^y ^^(^-e,(k))xnf[e,(k)]-^cos(z9)^. (14) 



Thus, as it is seen from (|T^), the ferromagnetie transition is linearly suppressed by the superexchange interaction, 
in agreement with previous results obtained iiJlil. The energy scale is defined in the units of the hopping amplitude. 
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I A \. The experimental results for the spin-wave stiffness coefficient provide | A |=0.16 eV (see e. and references 
therein) . 

As we see in Fig. 2, Tc decreases as system is doped with holes. For concentrations, where the CMR effect is 
very pronounced {x ~ 0.3), the value we found for the critical temperature is equal to 390 K, which is in good 
agreement with the experimental data. Hence, the DE model via degenerate orbitals alone gives a good value for the 
critical temperature. For small concentrations our result does not agree with experimental observations, in which Tc 
is iiicreasing with doping. Such a behavior for small concentrations can be explained in terms of percolation theory 
(seeErO and references therein) and can not be captured by using the present approach. 
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FIG. 2. Curie temperature is plotted as a function of hole concentration x. 

In conclusion, we have shown, that the variational mean field theory can successfully describe the paramagnetic- 
to-ferromagnetic transition for x ~ 0.3 in the frame of the DE model when the two-fold degeneracy of the eg levels 
is taken into account. We have also discussed the competition between the ferromagnetic and antiferromagnetic 
exchange interactions. Although the phonons might play an essential role in the high-temperature properties of the 
manganites near the metallic region (isotope effect), the other mechanism for suppressing the ferromagnetic phase is 
the superexchange interaction between the t2g core spins. 

We thank L. P. Gor'kov for very useful discussions and F. Drymiotis for help with the manuscript. This work was 
supported by the National High Magnetic Field Laboratory through NSF Cooperative Agreement ^ DMR-9527035 
and the State of Florida. 
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